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ABSTRACT 

We present a study of a sample of LMC red giants exhibiting Long Secondary 
Periods (LSPs). We use radial velocities obtained from VLT spectral observations and 
MACHO and OGLE light curves to examine properties of the stars and to evaluate 
models fo r the cause of LSPs. This sample is much larger than the combined previous 
studies of iHinkle et~aTl (|2002h and Wood, Olivier & Kawaler (2004). 

Binary and pulsation models have enjoyed much support in recent years. Assuming 
stellar pulsation, we calculate from the velocity curves that the typical fractional 
radius change over an LSP cycle is greater than 30 per cent. This should lead to 
large changes in T Q g that are not observed. Also, the small light amplitude of these 
stars seems inconsistent with the radius amplitude. We conclude that pulsation is not 
a likely explanation for the LSPs. The main alternative, physical movement of the 
star - binary motion - also has severe problems. If the velocity variations are due to 
binary motion, the distribution of the angle of periastron in our large sample of stars 
has a probability of 1.4 x 10 -3 that it comes from randomly aligned binary orbits. 
In addition, we calculate a typical companion mass of O.O9M0. Less than 1 per cent 
of low mass main sequence stars have companions near this mass (0.06 to 0.12 M©) 
whereas ~ 25 to 50 per cent of low mass red giants end up with LSPs. We are unable 
to find a suitable model for the LSPs and conclude by listing their known properties. 
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1 INTRODUCTION of these stars is usually found on sequence B (Wood et al. 

1999). 

A subset of Long-Period Variable stars (LPVs) were found 

several decades ago to show a secondary period of varia- Approximately 25-50 per cent of Long-Period Variables 

tion, in addition to their primary pulsation (e.g. ESS show an LSP (Wood et al. 1999; Sos z yriski et al. 2004b ; 

|Pavne-Ga D oschkirJ[l95lh . These Long Secondary Periods, | Percy et al. |[200l [ Soszynski et al. | » | Fraser et al. ||2008|). 

or LSPs, exceed the primary period in length by approx- At P resent > there is no accepted explanation for LSPs. Wood 

imately one order of magnitude, a fact more recently con- et aL ( 1999 ) initiall y Prosed a model wherein the sequence 

firmed by Wood et al. (1999). Wood et al. showed that LPVs D stars are semi-detached red giant binaries, as the length of 

fall on four distinct period-luminosity sequences A-D, and their P eriods are consistent with those expected for binary 

they found an additional sequence, E, of red giant binaries. s y stems with solar-mass components. Several other models 

These multiple sequences have been confirmed in subsequent have been P r °P° s ed to explain LSPs, includin g radial and 

studies, and a splitting of sequence B into tw o sequence s has nmiradial pulsation, and dust effects. However, | Wood et al. | 

since been di scover ed (Soszyriski et al. 2004a; llta et alJliool : » demonstrate that there are problems with all of these 

iFraser et allB : ISoszvnski et alJB . Long Secondary models - Tt is clear that further investigation is required in 

Periods occupy sequence D, the sequence corresponding to order to dlscover the cause of LSPs ' 

variations with the longest period. The primary pulsation Here > we stud y a sample of sequence D stars in the Large 

Magellanic Cloud (LMC) in order to provide new constraints 
on models for the sequence D phenomenon. In particular, 

* E-mail: nicholls@mso.anu.edu.au (CPN); we present new spectra taken with the ESO's Very Large 

wood@mso.anu.edu.au (PRW); m.cioni@herts.ac.uk (M-RLC); Telescope (VLT). We derive velocity curves, radii and effec- 

soszynsk®. astrouw.edu.pl (IS) tive temperatures from these spectra and compare them in 



© 2009 RAS 



2 C. P. Nicholls et al. 



phase with light curves and associated quantities from the 
MACHO and OGLE databases. 

Our study f o llows the methods of iHinkle et al.l (|2002l ) 
and IWood et all l|2004) . who also compared both light and 
radial velocity data for sequence D stars. However we use a 
larger dataset and have the advantage of high quality and 
simultaneous light and velocity data, which we hope will 
give a more complete and accurate picture of the behaviour 
of LSP variables. Since our sample comes from the LMC, 
we also have the advantage that the distances, and hence 
derived luminosities, are well-determined. 



2 OBSERVATIONS AND DATA REDUCTION 

We chose a region of high star density in the middle of 
the LMC bar for this study, in order to get a large sam- 
ple of sequence D stars in the 20 arcmin field of the 
FLAMES multifibre system on the VLT. The field centre 
is at 05' 1 28'15" - 69°45'43" J2000. It contained a sample 
of 58 variable red giants exhibiting LSPs which were ob- 
tained from the MACHO Project database. It was impor- 
tant to have simultaneous light and velocity data. Fortu- 
nately, OGLE light curves were taken at the times of our 
VLT observations. We have used both OGLE II and OGLE 
III I-band light curves in this study. 

Spectra were obtained on 21 nights from 2003 Novem- 
ber to 2006 March. The FLAMES/ GIRAFFE spectrograph 
(Pasquini et al. 2002) with a grating setting of HR16 was 
used, giving spectra with a resolution of 23900 and a spec- 
tral interval of 693.7-725.0 nm. This region includes the TiO 
bandhead at 705 nm, which can be used for spectral typing 
and to examine variation of T e g, as well as derivation of 
radial velocities. An exposure time of 20 minutes was used 
for all spectra. The spectra were obtained in service mode. 
Unfortunately no observations were taken between May and 
August when the Large Magellanic Cloud (RA = 05 h 23') 
was only observable near twilight, leaving a long Winter gap 
in 2004 and 2005. 

The raw spectra of the 58 sequence D stars, most of 
which have 21 observations from different dates were re- 
duced using the FLAMES/GIRAFFE pipeline, or the Swiss 
reduction pipeline. A number of sample spectra are shown 
in Fig. [TJ where the TiO bandhead is clearly visible in the 
second and third panels. The majority of the program stars 
have spectra similar to those shown in the top three panels, 
i.e., they are oxygen-rich giants of spectral types M and K. 
However, the sample does contain a number of carbon-rich 
stars (C-stars), as shown in the fourth panel of Fig. [TJ 

Using the iraf software, each spectrum was visually ex- 
amined and any obvious cosmic rays removed. There was a 
small percentage of bad data: in particular, spectra taken 
on two dates (2005 March 17 and 2005 September 10) had 
very few counts. This was due to the weather and seeing 
conditions at the time, and the spectra taken on these dates 
were discarded, leaving 19 good observations for most stars 
(some stars had fewer spectra due to the fibre allocation 
process in FLAMES). A small number of other spectra were 
discarded due to a low number of counts. These spectra did 
not fall on particular dates or particular target stars, and 
were distributed randomly throughout the dataset. Presum- 
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Figure 1. Sample spectra from the VLT. Top panel: The spec- 
trum of an O-rich program star with no TiO bands. Second panel: 
A program star with weak TiO bands. Third panel: A program 
star with strong TiO bands. Fourth panel: A C-rich program star. 
Bottom panel: A telluric star. Program stars are identified by 
their MACHO numbers. Note the different scale of the y-axis in 
the lower two panels. 



ably, the low counts were due to the fibres for these stars 
being misaligned. 

Relative radial velocities were obtained through cross- 
correlation with the iraf task fxcor. From each star, a spec- 
trum with a high number of counts and narrow lines was 
selected, and this acted as the template spectrum for the 
star's cross-correlation. The cross-correlation was performed 
in the wavelength region 6950-7160 A, as this region is rela- 
tively free of telluric lines, as can be seen in the lower panel 
of Fig. [T] Errors in the radial velocities were also taken from 
fxcor. The mean velocity error of the sample is 0.42km s _1 , 
and errors typically range between 0.1 and 0.6kms _1 . 

Because of uncertainty about the accuracy of the ve- 
locity calibration resulting from the pipeline reduction, the 
spectra were checked to see if any zero-point correction to 
the velocity was needed, before we calculated their absolute 
velocities. This procedure is now described. 



2.1 Testing the Relative Telluric Radial Velocities 

To check the velocity calibration, program stars were cross- 
correlated with a B star whose spectrum contained only tel- 
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luric lines. Cross-correlation of a program spectrum to a tel- 
luric spectrum produces the relative telluric velocity of the 
two spectra, which should be zero if there are no calibration 
errors. This cross-correlation was performed in the wave- 
length region 7160-7220 A, as it has a large concentration 
of telluric lines (see Fig. [TJ . 

For most dates of observation, the relative telluric ra- 
dial velocities were found to be not significantly different 
from zero, and so no systematic correction was applied to 
the spectra for those nights. However observations from two 
nights, with Heliocentric Julian Dates (HJD) of 2453377 and 
2453643, consistently showed a positive velocity offset from 
zero, so a zero-point correction for the velocities on these 
dates was included for all stars. The velocity corrections for 
these dates were — 2.061km s" 1 and -0.729km s" 1 respec- 
tively. 



2.2 The Effect of Telluric Lines on Relative 
Radial Velocities 

The highest concentration of telluric lines was found in the 
wavelength region 7160-7220 A, but there were a small num- 
ber in the program object region (6950-7160 A, see Fig. [T|. 
There was a possibility that these lines would affect the 
calculated radial velocities, so this was examined. The tel- 
luric lines in the program object region of several program 
stars were removed using iraf's telluric command. Then 
the program spectra with telluric lines removed were cross- 
correlated against each other in the same way the original 
program spectra were. It was found that removing the tel- 
luric lines in this region had no significant effect on the ve- 
locities or errors, so telluric line removal was not performed. 

2.3 Calculating the Absolute Radial Velocity 

Following the zero-point corrections, the absolute radial ve- 
locities of the program stars were calculated. First, the tem- 
plate of each program star was cross-correlated to a star of 
known radial velocity, giving the absolute radial velocity of 
each template. Then by adding a keyword 'vhelio' (which 
contained the absolute radial velocity) to the image header 
of the template, we could cross-correlate the spectra for each 
star with its template to obtain the absolute radial veloci- 
ties. 

For the O-rich program stars, the radial velocity ref- 
erence used was a Cet. Its spectrum was taken using the 
echelle spectrograph (resolution 70000) on the late 74-inch 
telescope at Mount Stromlo Observatory, Canberra, Aus- 
tralia. For the C-stars, the reference used was the C-star 
X Vel, for which a spectrum was also taken with the 74- 
inch. The radial velocity used for a Cet was — 25.8kms _1 
(taken from the Astronomical Almanac). The velocity for 
the spectrum of X Vel (— 5.4kms ) was obtained by cross- 
correlation with a Cet (even though the spectral types are 
different, there was a strong cross-correlation peak due to 
common metal lines). 

Radial velocities as a function of heliocentric Julian date 
(HJD) for a sample of the sequence D stars are shown in Ta- 
bleQ] Stars are identified by their MACHO numbers. Dashes 
denote dates for which there is no spectrum. The full table 
can be viewed in the online version of this paper. 



2.4 Correcting for the Short-Period Variation 

Because sequence D stars show two distinct modes of vari- 
ability, it is necessary to correct the light and velocity data 
for the behaviour of the short (primary) period in order 
to accurately study the LSP. We need simultaneous veloc- 
ity and light data for this purpose, so we use the OGLE I 
light curve in these procedures, even though it is missing the 
Southern Winter season. 

First the OGLE I light curve was interpolated and 
boxcar-smoothed over a time interval equal to the (primary) 
pulsation period. Once this LSP-only light curve was made, 
it was subtracted from the original light curve, leaving only 
the variations due to pulsation, di, or the 'pulsation light 
curve'. 

Similarly to obtaining the pulsation light curve, we ob- 
tained the pulsation velocity curve by making a binary fit 
to the velocity data using the period of the LSP. The devia- 
tions of the velocity data, dv, from this fit form the pulsation 
velocity curve. We emphasise that although a binary fit was 
used for the purpose of obtaining a smooth fit, binarity may 
not be the physical mechanism behind the LSP. 

For radial pulsation the light and velocity variations are 
shifted in phase relative to one another, so it is necessary to 
find the phase shift between the pulsation light and velocity 
curves in order to find the amplitude relation between the 
short-period light and velocity. 

The phase correction was found by plotting the pul- 
sation light versus pulsation velocity. By varying the phase 
shift of this plot, and judging by eye at which phase shift the 
pulsation light and velocity showed the clearest correlation, 
the best phase shift was found to be ~0.25 of a cycle. This 
means that most often the I light minimum occurs close to 
the same phase as the mean rising radial velocit y . Lebz elter, 
Kiss & Hinkle (2000) and ILebzelter fc Hinklei (|2002h find 
that the phase shift between the light and velocity varia- 
tions of pulsating SRVs is 0.5. This means maximum radial 
velocity occurs at minimum light, and vice versa. The phase 
shift is discussed further below. 

Next, the slope of the light-velocity amplitude correla- 
tion for pulsation was measured, to give the velocity correc- 
tion. From examination of several stars with relatively large 
pulsation amplitudes, it was concluded that the velocity cor- 
rection was ~15 km s _1 mag -1 in I. Generally, the correction 
to the LSP velocity resulting from these processes is small, 
less than 0.6km s _1 , as the median value of di is ~0.04 mag. 

In order to further investigate the phase relation be- 
tween the velocity and light curves, we constructed a model 
semi-regul ar pulsator usin g the n onlinear pulsation code de- 
scribed in lKeller fc Wood ( 20061 ). The model is a first over- 
tone pulsator, appropriate for a sequence B variable. The 
model oscillations should be similar to the primary oscilla- 
tion in the sequence D vari ables studied here, and the semi- 
regular variables studied by ILebzelter fc Hinkkj |2002). The 
model had M = 1.5M Q , L = 3000L Q , T eS = 3740^, helium 
abundance Y — 0.3 and metallicity Z = 0.004. My and 
Mi were computed from the M star model atmospheres of 
iHoudashelt et all (2000), assuming [Fe/H] — —0.5. 

The model light and radial velocity curves are shown in 
Fig. f2] The graph also shows the time variation of T e g and 
radius for comparison with observational estimates of these 
quantities in later sections. 
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Table 1. Heliocentric Radial Velocities of Sequence D Stars (kms 1 ) 



HJD 


77.7427.41 


77.7427.72 


77.7428.120 


77.7428.65 


77.7429.100 


77.7429.47 


77.7429.64 


77.7430.46 


77.7430.51 


2954.85 


240.85 


273.53 


273.61 


297.74 


314.55 


280.43 


278.15 


244.02 


266.31 


3005.86 


239.86 


274.41 


272.44 


299.37 


313.72 


280.73 


277.68 


245.33 


265.78 


3067.60 


241.32 


275.08 


272.69 


298.12 


313.20 


279.58 


277.01 


244.83 


267.10 


3091.56 


239.64 


271.16 


271.76 


298.36 


314.18 


279.08 


276.49 


244.71 


266.76 


3280.86 


241.59 


273.70 


272.70 


295.91 


314.24 


278.73 


- 


241.53 


264.68 


3324.76 


239.26 


275.30 


273.26 


294.60 


314.51 


278.14 


- 


240.90 


263.61 


3344.80 


240.04 


274.01 


273.34 


295.59 


315.15 


277.93 


- 


242.23 


263.88 


3376.58 


239.38 


274.15 


272.60 


295.73 


314.17 


277.91 


- 


243.88 


263.67 


3418.69 


238.11 


273.68 


272.59 


296.53 


315.73 


277.24 


- 


244.06 


265.99 


3471.51 


237.08 


271.80 


272.23 


298.19 


314.98 


278.01 


- 


243.95 


266.41 


3623.88 




















3642.83 


240.15 


271.86 


273.68 


297.61 


313.64 


281.31 




242.22 


267.86 


3644.82 


239.89 


271.50 


273.67 


297.36 


313.43 


281.03 




241.70 


267.60 


3663.86 


240.49 


271.74 


273.28 


296.62 


312.03 


281.83 




241.78 


268.17 


3664.86 


240.62 


271.92 


273.53 


296.83 


311.85 


281.99 




241.45 


268.16 


3685.67 


240.85 


273.02 


272.99 


297.64 


313.59 


282.93 




240.64 


267.92 


3709.79 


240.50 


269.51 


273.40 


297.15 


314.57 


284.06 




242.48 


266.48 


3741.60 


242.04 


272.47 


272.40 


296.61 


313.84 


282.54 




243.07 


265.69 


3768.62 


241.31 


273.48 


274.19 


297.47 


313.53 


280.68 




242.40 


266.54 


3819.50 


242.35 


274.62 


273.97 


296.13 


313.88 


282.12 




245.04 


266.56 



The model shows that minimum light occurs between 
mean increasing radial velocity (as estim ated here) and max- 
imum radial velocity (as estimated by iLebzelter fc H inkle 
2002). The model also predicts a ratio of velocity to I am- 
plitude of 9.2kms _1 mag -1 , somewhat smaller than the esti- 
mate of 15km s" 1 mag -1 given here. Given the small size of 
the corrections to velocity that we have applied, the uncer- 
tainty in the relative phase and amplitude ratio of velocity 
and I amplitude will not significantly affect the resulting 
velocity curve of the sequence D variability. 



2.5 Obtaining LSP Parameters and Narrowing 
the Sample 

Following the corrections described above, we were now in 
a position to analyse the light and velocity variations as- 
sociated with the Long Secondary Period. A binary fit was 
made to the corrected velocity data of the LSP using a For- 
tran program, Fitall to obtain the parameters of the ve- 
locity curve. Fitall also made a Fourier series fit (with a 
frequency f and one harmonic 2f) to the MACHO Mb, Mr 
and OGLE / light curves, and to {Mb — Mr). Some exam- 
ples of the fits to the light and velocity data made by Fitall 
can be seen in Fig. [3] 

Due to the poor quality of the data for some stars, the 
original sample of 58 sequence D stars had to be reduced. 
Some stars had large velocity errors due to poor-quality or 
few spectra or no OGLE light curve, so they were excluded 
from the sample. A few stars had LSPs near 365 days leading 
to a large unfilled gap in the velocity curve, and could not 
be analysed. 

Next, for each of the remaining stars we judged by eye 
the quality of each star's binary fit, and eliminated all stars 
judged to have poor fits to their velocity data. Finally, we 
were left with a reduced sample of 30 sequence D stars, from 
the original 58. The parameters of the velocity fit for these 
30 stars are shown in Table[2] where 7 is the system velocity, 



K is the velocity semiamplitude, e is the eccentricity, ui is 
the angle of periastron, T is the date of periastron, P is the 
period, and f(m) is the mass function, given by 



f{m) = 



K P 



(1) 



2ttG (m + A/) 2 

in which M is the mass of the star, and m the mass of its 
companion. 



3 RESULTS 

In this section, we discuss some general results that are not 
model specific. Model specific results are presented in sec- 
tion m 



3.1 The True Distribution of the Velocity 
Amplitude 

The distribution of velocity amplitude is an important pa- 
rameter for most models of the sequence D phenomenon. 

A histogram of the velocity amplitudes of our sample 
is shown in Fig. [4] with the values separated into bins of 
width 0.5km s . It is important to note that we use the 
full or peak-to-peak amplitude, as opposed to the semi- 
amplitude. These amplitudes are obtained from the binary 
orbit fits to the observed velocities. We find that the ve- 
locity amplitude is concentrated mainly between 2.5 and 
5.0km s" 1 , with the highest number of stars falling in the 
bin centred at 3.5km s -1 . The distribution has a median 
value of 3.53km s _1 . Thus our sequence D stars all have 
low, relatively similar velocity amplitudes. These results 
are consistent with the results of Irlinkle et al.l ((2002) 's and 
IWood et al.l (|2004t )'s analysis of Galactic LSP variables. 

There are very few stars in our sample with velocity 
amplitudes < 3.0km s _1 . To find out if we could detect such 
stars we performed a Monte Carlo simulation. Fifty stars 
were simulated for each of the velocity amplitudes of 3.0, 2.0, 
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Figure 3. Light and velocity curves plotted against phase for four sequence D stars. The blue velocity points without error bars are the 
radial velocities without corrections and the red points with error bars are the radial velocities that have been corrected for pulsation 
and had specific date corrections. Also shown are the Fourier fits to the light data (black lines), and the Binary fits to the velocity data 
(cyan lines), both made by FlTALL. 
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Table 2. Orbital Elements 



Star 7 K e u TP f(m) No. 

(kms" 1 ) (kms" 1 ) (deg) (HJD) (days) (Kr 4 M Q ) obs. 



77. 


7427. 


.41 


239 


.812 


± 


0. 


.144 


1. 


.870 


± 


0. 


175 


0. 


436 


± 





.139 


189.68 ± 


6 


.91 


3491.19 ± 13 


.82 


776. 


18 


3 


.835 


± 


1 


.377 


19 


77. 


7427. 


.72 


273 


.087 


± 


0. 


.186 


1. 


.626 


± 


0. 


.271 


0. 


547 


± 





.182 


257.71 ± 


6 


.59 


3251.47 ± 10 


.27 


488. 


67 


1 


.277 


± 





.840 


19 


77. 


7428. 


.65 


297. 


.159 


± 


0. 


.116 


1. 


.850 


± 


0. 


095 


0. 


402 


± 


0. 


.115 


290.72 ± 4 


.53 


3462.77 ± 4 


.97 


498. 


94 


2 


.512 


± 


0. 


.567 


19 


77. 


7429. 


.100 


313 


.955 


± 





.145 





.966 


± 


0. 


.188 


0. 


187 


± 





.145 


182.16 ± 12 


.49 


3166.86 ± 16 


.46 


473. 


24 





.419 


± 





.248 


19 


77. 


7429. 


.47 


280. 


.326 


± 





.108 


2 


.271 


± 


0. 


.120 


0. 


.475 


± 





.108 


274.88 


± 3 


.59 


3547.93 ± 9 


.50 


1078. 


33 


8 


.910 


± 


2 


.264 


19 


77. 


7430. 


.46 


243. 


.159 


± 


0. 


.129 


1. 


.687 


± 


0. 


161 


0. 


195 


± 





.127 


279.99 ± 5 


.57 


3370.88 ± 6 


.56 


416. 


45 


1 


.957 


± 


0. 


.579 


19 


77. 


7430. 


51 


266. 


.274 


± 





.212 


2 


.203 


± 


0. 


137 


0. 


583 


± 





.202 


227.13 ± 3 


.50 


4148.37 ± 11 


.12 


798. 


34 


4 


.737 


± 


2. 


.686 


18 


77. 


7547. 


.41 


271. 


.204 


± 


0. 


.148 


2 


.150 


± 


0. 


163 


0. 


151 


± 





.147 


87.27 ± 4 


.59 


3314.19 ± 10 


.26 


870. 


06 
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1.0 and 0.5km s -1 . The period, eccentricity, and time of pe- 
riastron were all selected randomly. The number of velocity 
points in each simulation was the same as the number of ob- 
servations (19) and they were distributed in time identically 
to the observations. Random noise was added to each point, 
drawn from a normal distribution whose mean was the mean 
velocity error in the observations. In plots for visual exami- 
nation, an errorbar was added for each point, drawn from a 
uniform distribution between 0.2 and 0.5km s _1 (consistent 
with observed errors). 

A binary fit was made to the simulated velocities in the 
same way as to the observed velocities (see section [23}. The 
simulated data were put through the same reduction steps as 
the observed data, to see what percentage would pass. The 
percentages can be seen in Fig. [5] This tells us that we should 
be able to detect velocity amplitudes < 3.0km s _1 . Using our 
simulated percentages, we can estimate the true number of 
stars at each velocity amplitude by dividing the observed 
histogram by the observable percentage. The resulting mod- 
ified velocity amplitude histogram is shown overplotted on 
the observed histogram in Fig. [4] 

The fact that we do not observe many sequence D stars 
with velocity amplitudes < 3.0km s _1 suggests there are not 
many to be observed. We conclude that the true distribution 
of sequence D velocity amplitudes, represented in Fig. [4] 
shows a strong peak near 3.5km s . 



3.2 Correlation of Light and Velocity Amplitudes 

Light amplitude is plotted against velocity amplitude in 
Fig. [5] There is no obvious correlation between light and 
velocity amplitudes in either Mb, Mr, or I. As noted previ- 
ously the velocity amplitude has a small spread, and peaks 
at the small value of 3.5km s . In particular, stars with 
small light amplitudes have velocity amplitudes as large as 
those stars with the largest light amplitudes. 



3.3 Correlation of Velocity Amplitude with 
Magnitude and Period 

Velocity amplitude is plotted against K magnitude in Fig. [7] 
which shows a suggestion that velocity amplitude increases 
slowly with luminosity. A similar correlation is seen in Fig. [8] 
between velocity amplitude and LSP. However these possible 
correlations are disrupted by the presence of a few stars with 
larger velocity amplitudes, that fall apart from the general 
relation. 

4 TESTING AND CONSTRAINING MODELS 

Here we examine some of the more plausible and popular 
models of LSPs, in the light of our new data. 
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Figure 2. Mj, ;, Mj, My, R/Rq, T e g and v r plotted against 
time in a model red giant variable pulsating in the first overtone 
mode. The radius R is defined at optical depth ^ and v T is the 
radial velocity seen by a distant observer. It is defined as —v/lA, 
where v is the pulsation velocity, relative to the centre of the star, 
of matter near optical depth g. The factor 1.4 is the correction 
factor from observed t o pulsation velocities for red giant pulsators 
l|Scholz fc Woodll2000h . 



Figure 4. A histogram of the observed full velocity amplitude 
is plotted in red (forward shading). Note the clustering around 
S.Skms^ 1 . Overplotted in blue is the modified histogram to ac- 
count for the observational detection probability calculated in our 
Monte Carlo simulation (see text for details). 



4.1 The Binary Model 

The model of binary motion as an explanation for LSPs 
has been much debated since its proposal by Wood et al. 
(1999). In its favour, binary motion may be able to explain 
the sequence D velocity curves, and the light variation was 
originally thought by Wood et al. to resemble an eclipse 
of the star by a dusty cloud surrounding a small, orbiting 
co mpanion. How e ver w hen t his model wa s furth er examined 
by iHinkle et all (|2002h and IWood et all (|2004l ). two main 
problems were found: the companion mass was calculated to 
be always ~0.1 Mq, and the characteristic shape of sequence 
D stars' velocity curves was shown to imply that the angle 
of periastron is not uniformly distributed between and 2-k. 
With our larger sample we can place tighter constraints on 
these properties. 



4-1.1 The Companion Mass 

In a binary system, the velocity amplitude and period give 
an estimate of the mass of the companion. For a typi- 



cal sequence D star with an LSP of 500 days, a typi- 
cal velocity amplitude of 3.5 kms -1 , and an assumed to- 
tal system mass of M = 1.5 Mq, the orbital separa- 
tion is of the order of 1.4 AU and the companion has 
a mass of 0.09 Mq. This means it is a very low mass 
main sequence star or a brown dwarf. There is an ob- 
served deficit of brown dwarf companions to main sequence 
stars, compared to both more massive stellar companions 
and less massive planetary companions. This is known as 
the 'Brown Dwarf Desert' llMcCarthv fc Zuckermanl 120041 : 
iGrether fc Lineweaverl [2006 ) . The range of likely compan- 
ion masses can be found by calculating the companion 
mass for velocity amplitudes within one sigma of the mean 
of 3.5km s . We find a companion mass range of 0.06 
to 0.12 Mq. Using the fit s to the data in Fig. 8 of 
IGrether fc Lineweaverl (12006! ). we find that only 0.86 per 
cent of low mass main sequence stars have companions with 
masses in this range. Thus binarity is an unlikely model, 
when we remember that ~30 per cent of low mass stars will 
exhibit LSPs when they pass through the AGB stage. 
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Figure 5. Fraction of simulated velocity curves which pass our 
data reduction process, for different velocity amplitudes. 



Figure 6. Light amplitude plotted against velocity amplitude for 
M B , Mr, and /. 



4-1.2 Correlation of Velocity Amplitude with Light 
Amplitude, Magnitude and Period 

We would not expect the velocity and light amplitudes to 
be related in a binary model where the light variation is 
due to an eclipse phenomenon (by a cloud in the Wood et 
al. 1999 model). However, these amplitudes may be related 
if the light variations were caused by distortion of the red 
giant by its unseen companion (see section I4.1.5P . 

As noted in section 13.31 only a hint of a correlation is 
seen between velocity amplitude and K magnitude (Fig. [7}, 
and with LSP (Fig. |SJ). We cannot draw any conclusions 
from this in the context of a binary model. 



4-1.3 The Distribution of the Angle of Periastron 

One of the most conclusive ways to test the plausibility of 
the binary model as a cause for the LSPs is to examine the 
distrib utio n of the angl e of p eriastron, u>. As iHinkle et al.l 
l|2002h and IWood et all (|2004 ) note, the observed distribu- 
tion of u> should be uniform over the whole range of angles, 
since one would expect binary orbits to be randomly aligned 
in space. 

The angle of periastron was calculated by the program 
Fitall as part of the binary fit made to our radial velocity 
data. The distribution of the angle of periastron is plotted 
as a histogram in Fig. [9] with the data divided into bins 



20° in width. Note that since the errors in the angle of pe- 
riastron are typically around 5° (Table [2}, the distribution 
of angle of periastron shown in Fig. [5] will not be signifi- 
cantly broadened by these errors. There is a clear bias in 
the distribution towards angles > 180°, and the median of 
the sample is 227°. The distribution of the a ngles of pe- 
riastro n calculated bv lHinkle et al. (2002) and lWood et al.l 
(2004) for their small samples of Galactic stars has been 
overplotted. It is clear that their values of uj lie in the same 
region as the bulk of the present sample. This shows that 
the LSPs share a common cause, whether found in the LMC 

or i n our own Ga l axy. 

IHinkle et all l|2002ft and I Wood etU] (|2004h both found 
that the characteristic shape of the velocity curves suggests 
an eccentric orbit and a large angle of periastron, and that 
the distribution of the calculated angle of periastron was in- 
consistent with what would be expected if the sequence D 
stars were binaries, as no satisfactory explanation could be 
found for this bias towards higher angles. But given their 
small sample size, their conclusions were not highly signifi- 
cant. We can quantify and greatly improve these claims by 
using a Kolmogorov-Smirnov, or K-S Test, with our new 
data. 

We have used the one-sample K-S Test, in the form of 
the F ortran 77 subroutine ksone (provided in iPress et al.l 
1986), to compare our distribution of the angle of periastron 
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Figure 7. Velocity amplitude plotted against K magnitude. 



Figure 8. Velocity amplitude plotted against LSP. 



to the uniform distribution. We find that the probability 
that our distribution of u> is consistent with the uniform 
distribution is 1.4 x 10 . In other words, the probability 
that the sequence D stars are binaries is extremely small. 
This is a major result from this study. 

An angle of periastron between 180° and 360° means 
that at periastron, the red giant is closest to the observer, 
with the smaller companion further away. It is hard to see 
what selection effects could cause an LSP to be observed 
only in these binary orientations. 



4-1-4 The Distribution of the Eccentricity 

A histogram showing the distribution of the eccentricity is 
shown in Fig. 1101 The eccentricity was calculated by Fitall 
as one of the orbital parameters of the binary fit to the 
velocity data. 

The plot shows that if the sequence D stars are caused 
by binary motion, then they are in eccentri c orbi ts. This con- 
firms the sug gestion of iHinkle etafl (|2002u and I Wood et al.l 
(|2004l V The median eccentricity of the sample is 0.3. 



4-1.5 Ellipsoidal Variables 

An ellipsoidal variable is a binary star which is tidally dis- 
torted by an orbiting companion into an ellipsoidal shape. 
Light variability is caused mostly by the variation of the 



apparent surface area of the star seen by the observer as 
the star rotates. Its velocity curve is dominated by orbital 
motion, with a small contribution from the rotation of the 
limb-darkened ellipsoid. 

The light curve for an ellipsoidal variable is expected to 
have two maxima and minima per orbit, due to the star's 
shape. However, the velocity curve is expected to have only 
one maximum and minimum per orbit. We can use this char- 
acteristic to investigate the plausibility of ellipsoidal vari- 
ability as an explanation for the LSP variation. 

Soszyhski et al. (2004b) show that the light curves of se- 
quence E variables are satisfactorily explained by ellipsoidal 
variability. Their data show a partial overlap of sequences 
D and E in the P-L diagram, and so they suggest that the 
Long Secondary Periods may have a binary ellipsoidal origin. 
ISoszvhskil (2007) explores this possibility further by search- 
ing for ellipsoidal or eclipsing shapes in residual LSP light 
curves. He finds these variations in ~5 per cent of his sam- 
ple, and adopts the binary model to explain LSPs. 

However, with our new dataset (which includes some 
sequence E stars), the difference between sequence D and 
sequence E mechanisms can now be vividly demonstrated 
by the difference between the light and velocity curves of 
sequence D and E stars. The phased-up light and velocity 
curves of sequence E stars show that the light curve com- 
pletes two full cycles to every single cycle of the radial ve- 
locity, as expected for an ellipsoidal variable (|Adams et al.l 
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Figure 9. The distribution of the angle of periastron, u). This 
study's sample i s plotted in r ed (f orwar d shading). The co m- 
bined sample of Iffinkle et all ||2002| ) and IWood et all (|2004|) is 
shown in blu e (back shading). Two s tars have bee n dropped from 
Hiiikl e~et all s sample and one from lWood et all 's, due to overly 
large errors for w. 

2006, Nicholls, Wood & Cioni, in preparation). However, 
Fig. [3] shows that sequence D stars do not display this be- 
haviour: the phased-up light and velocity curves match each 
other cycle for cycle. Furthermore, sequence E stars typi- 
cally have velocity amplitudes of 30km s" 1 or more while 
the sequence D st ars have significan tly lower typical ampli- 
tudes of 3 . 5km s ~ 1 . ISoszvi iski (2007) predicted that sequence 
D stars showing residual ellipsoidal or eclipsing-type varia- 
tions should have large velocity amplitudes, similar to those 
observed for s equence E stars. One of our stars (77.7671.282) 
was noted by ISoszvriskil as a double-humped LSPV and its 
velocity amplitude is only ~4 kms -1 . Additionally, none 
of our LSP sample show typical sequence E velocity ampli- 
tudes (see Fig. 2} . Sequence D and E variables clearly have 
distinctly different velocity amplitudes. 

There is an unambiguous distinction between the mech- 
anisms of variability responsible for causing sequences D and 
E. We can now state with some certainty that binary star 
ellipsoidal variability is not the cause of the LSP. 

Similar to the sequence E ellipsoidal variables are sym- 
biotic binaries, which are red giants with an accreting white 
dwarf companion, often in an eccentric orbit. Radial veloc- 
ity amplitudes of these systems are usually ^ 8km s _1 , and 



12 



10 - 




0.2 0.4 0.6 0.8 

Eccentricity 



Figure 10. The distribution of the eccentricity. 

the red giant may be a Mira variable l|Hinkle et al.l I2006T ) . 
However, symbiotic spectra show high temperature emission 
lines from accretion onto the white dwarf, something that 
has not been seen in either the sequence E or sequence D 
spectra. There does not seem to be any relation between the 
symbiotic stars and the sequence D stars. 

4.2 Radial Pulsation 

There are several ways in which the present dataset can be 
used to test radial pulsation models: 

4-2.1 The Phase Difference between Light and Velocity 
Curves 

Minimum light and the mean of the velocity during increase 
are relatively well defined for our sequence D stars. In order 
to examine their relative phases, a histogram of the phase 
difference between minimum light and mean rising velocity 
is shown in Fig. 1111 The phases of minimum light and mean 
rising velocity were obtained from the fits made to the light 
and velocity curves. Each panel shows a strong peak of phase 
differences lying between and 0.2 and there is little spread 

in p hase. 

iLebzelter efaLI (2000) and [Ceb zelter fc Hinklel (2002) 
find that for radially pulsating SRVs, the light lags the ra- 
dial velocity by roughly 0.5 in phase, which means that min- 
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Figure 11. Histograms of the difference between the phase of 
LSP minimum light and the phase of mean rising LSP radial 
velocity. The top panel shows the phase difference for MACHO 
Mb, the middle panel Mr, and the bottom panel OGLE /. Radial 
velocity is defined with respect to the observer: a positive radial 
velocity corresponds to matter moving away from the observer 
i.e. infalling to the centre of the star in a radial pulsation model. 

imum light occurs ~0.25 in phase after the time of mean ris- 
ing velocity. Our model SRV (Fig. shows that minimum 
light occurs close to 0.23 in phase after mean rising veloc- 
ity. Our sequence D stars are significantly different. Taking 
the main peak of light-velocity phase differences in Fig. 1111 
(values between —0.1 and 0.3), the mean phase shift is 0.1, 
and the standard deviation is 0.08. More simply, minimum 
light usually occurs ~0.1 in phase after mean rising velocity. 
Radial pulsation therefore may not be consistent with this 
aspect of the LSP phenomenon, although this is not a strong 
conclusion. 

4-2.2 The Shape of the Velocity Curve 

During large amplitude radial stellar pulsation (as in Mira 
variables) , a shock wave from the interior of the star causes 
the surface layers to rapidly accelerate outwards, and then 
settle back more slowly. This is also shown in the shape of 
the velocity curve of our low amplitude model SRV, Fig. [2] 
The velocity curves of most of our sequence D stars show ex- 
actly the opposite of this behaviour: typically we find a rapid 
increase in radial velocity with a correspondingly slower de- 



cline (see Fig. [3} . In a radial pulsation model, this translates 
to a slow increase in stellar radius and a quick decrease. This 
has not been observed in any r ed giant stars known to be ra- 
dially pulsating variables (e.g. lLebzelterlll999l ; iHinkle et al.l 
1 1984 for SR variables and Mira variables respectively). 

4-2.3 Correlation of Velocity Amplitude with Light 
Amplitude, Magnitude and Period 

The lack of correlation between velocity amplitude and light 
amplitude for the stars in our sample (see Fig. [6} is an un- 
expected result for pulsation mechanisms as they would be 
expected to show a positive correlatio n between light and 
velocity amplitude (|Hinkle et al.|[l997h . Additionally, corre- 
lations of velocity amplitude with luminosity and with pe- 
riod are expected for radially pulsating variables. However 
Figs.[7]and|8]show that any correlations between these quan- 
tities in our sample are vague at best. We are unable to draw 
any conclusions from this in the context of radial pulsation. 

4-2.4 The Variation of Stellar Radius 

The stellar radius can be calculated from our data in sev- 
eral ways. Our first method was to calculate R using the 
Stefan-Boltzmann equation L = A-kR 2 aT^ s , with L and T e g 
calculated from MACHO photometry (see Appendix). 

Examples of the variation of stellar radius calculated 
from the photometry, i? p hot , can be seen in Fig. 1121 There 
is generally a clear periodic curve at the period of the LSP 
as well as a scatter due to the primary pulsation. 

Variation of stellar radius can also be calculated from 
the velocity, using AR = J vdt where v is the fit to the ve- 
locity data made by Fitall. If we assume that the velocity 
changes are due to radial pulsation, then this gives the ex- 
pected radius variation for radial pulsation, R ve i, where the 
resulting radius has been normalised to the median photo- 
metric radius. Here the velocity fit data was scaled by a fac- 
tor of 1.4 to c onvert from observed to true pulsation ve locity 
in red giants l|Scholz fc Woodll2000l : IWood et ai]|2004h . The 
radius change computed from the velocity fit is shown in 
Fig. 1121 Again it shows a clear periodic curve. 

It is also possible to calculate stellar radius from the 
Stefan-Boltzmann equation using an effective temperature 
derived from the spectra, as opposed to the photometry 
(spectral T e « derivation is described in section I4.2.5|l . The 
luminosities at the times of the spectra were in this case 
calculated from the fits to Mb, Mr and (Mb — Mr) as 
described in the Appendix. The radius calculated in this 
manner, i? S p CC , can also be seen in Fig. 1121 

A histogram comparing the median values of -Rphot and 
fispec is shown in Fig. 1131 For both radius estimates, most 
stars have median radii lying between 100 and 200 R©. The 
median values for the entire sample for i? p hot and Rs pcc are 
135.4 and 132. 1R©. _R ve i is not shown, as it is normalised to 

-Rphot- 

The radial amplitudes of our stars mostly lie between 3 
and 60 R© for all three variations. The median amplitudes 
of the sample for AJ? p hot, A7? spcc and Ai? vc i are 5.92, 5.62, 
and 40.02Rq respectively. 

Fig. Q3] shows the distribution of radial amplitude AR 
divided by period, P. This quantity is a more sensitive in- 
dicator of differences between stars and models than the 
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Figure 12. Different properties of four typical sequence D stars plotted against Julian Date. Solid symbols denote values plotted at the 
date of measurement while open symbols are values that have been shifted forward or backward by one or more periods. Top panel: 
MACHO Red light curve. Second panel: MACHO Blue light curve. Third panel: OGLE / light curve. Fourth panel: Variation of J? p hot 
(small black points), variation of R vc i (cyan curve) and variation of R BpC c (large magenta points). Fifth panel: Variation of T a ff( phot ) 
(small black points), variation of T c ff( apec ) (large magenta points) and variation of T ef j( ve i) (cyan curve). Bottom panel: Observed radial 
velocity. See text for details. Plots for the entire sample are available in the online version of this paper. 
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Figure 13. Distribution of median values of -Rphot (solid red line) 
and i? S pcc (dashed blue line). 



Figure 14. Distribution of radial amplitude/period for i? p hot 
(solid red line, forward shading), R Bp cc (dashed blue line) and 
R ve l (dot-dashed green line, back shading). 



amplitude itself. In a radially pulsating model, one would 
expect the amplitude, AR, to be tightly correlated to the 
period if the velocity amplitudes are the same i.e. AR/P 
would be a constant. Given that our stars show very similar 
velocity amplitudes (see Fig. 3} we should expect tight, and 
similar, peaks for all three radius variations. 

From Fig. [T3] it is clear that the distributions of 
A7? p hot/f and Ai? spcc /P, are similar, but they are not the 
same as for AR ve \/P. Ai? ve i is computed assuming radial 
pulsation. Given that it appears that the three radius varia- 
tions do not agree, the assumption that AT? is due to radial 
pulsation of the star is unlikely. 

Fig. 1 151 shows the distribution of AR divided by R. For 
Pphot and -Rspcc, the fractional radius change is usually 10 
per cent or less. For i? ve i, the fractional radius change is 
much larger, usually between 10 and 50 per cent. Once again, 
the assumption that the measured radial velocity variations 
are due to radial pulsation does not seem to be correct. 

A histogram comparing the phases of the three radius 
variations to the phase of the Mr light curve is shown in 
Fig. 1161 We compared the phase of minimum radius with the 
phase of minimum light, as these are both well-defined parts 
of the curves. In order to calculate these phases a Fourier se- 
ries with four terms was fit to the radius curves. For the light 
curve, the previously-used Fourier fit was used (see Fig. [3]). 
Fig.[T6lshows that -Rphot and i? ap oc both vary almost in phase 



with the light variation (phase differences generally ^0.1). 
On the other hand, i? ve i is more likely to be out of phase 
with the light (and consequently with the other radius vari- 
ation estimates). Again, the assumption that the observed 
radius variations can be attributed to radial pulsation leads 
to inconsistencies. 

Using our model SRV (see section l2~4"ll minimum radius 
occurs about 0.7 in phase before minimum light (Fig. [2}. Of 
our variations, R ve i is the most similar to this. -Rphot and 
Rspec, whose phases agree well with each other, do not show 
the radius-light phase relation expected for radial pulsation. 

Overall, the discrepancy between the variations of R vc i 
and the other two radius estimates suggests that radial pul- 
sation is unlikely to be the cause of the LSPs. 



4-2.5 The Variation of Effective Temperature 

The variation of T e ff was calculated from the photometry, 
from the spectra and from the velocity fit. 

The variation in the photometric T e g (T e ff(phot) ) > whose 
derivation is described in the appendix, is shown in Fig. 1121 

T fr was also calculated using the depth of the TiO 
bandhead at 7054 A (T e ff( spec )). The depth was calculated 
as the ratio of the mean count level of a 4 A band long- 
ward of the bandhead to the mean count level of a 4 A band 
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Figure 15. Distribution of relative radial amplitude for -Rphot 
(solid red line, forward shading), -R S pec (dashed blue line) and 
-R ve l (dot-dashed green line, back shading). 



Figure 16. Distribution of phase of minimum light - phase of 
minimum radius for R p hot (solid red line, forward shading), R Bpe c 
(dashed blue line) and R VE i (dot-dashed green line, back shading). 



shortward of the bandhead, taking into account the fact that 
the LMC has an average velocity of 270km s _1 , leading to a 
longwards wavelength shift of 6.3 A in the spectra. 

In order to calculate the effective temperature from the 
TiO depth we calculated the depth of the TiO bandhead for 
9 Wing Spectral Standard stars. Then a quadratic fit was 
made to the TiO Depth-Spectral Type data. This quadratic 
fit was used to calculate the spectral types of our sequence D 
stars from their TiO depths, an d the effective t e mper atures 
were then calculated using the iRidgwav et all (|l980h 
Spectral type relation. An example of the typical variation 
of Tefffspoc) in these stars can be seen in Fig. 1121 

Some of the warmer sequence D stars in our sample 
showed little to no evidence of TiO in their spectra, and so 
we were unable to include these in the T e ff(s P oc) calculation. 
However, good-quality velocity curves were not necessary 
for calculating T cf f( spcc ) so a number of stars previously dis- 
counted from the velocity studies were added in where con- 
sideration of radius and effective temperature was required. 
Overall there was a net increase in sample size from 30 to 
37 (the stars added in were excluded only from the veloc- 
ity amplitude, angle of periastron, eccentricity and velocity 
phase parts of this work). 

We also calculated the effective temperature variation 
expected for radial pulsation, T eS ^ vel y This was calculated 



from the radius change computed for pulsation in sec- 
tion l4,2.4l the luminosity calculated from the fits to the pho- 
tometry, and the Stefan-Boltzmann equation. The variation 
of T e ff( vo i) can be seen in Fig. [TJ] 

The distributions of the median values for each star 
of T e fF( phot ) and T e ff(s P ec) are shown in Fig. [TTJ Most stars 
have median effective temperatures lying between 3300 and 
3800-K\ The median values for the entire sample for T e ff( P hot) 
and T c ff( spcc ) are 3662.9 and 3725. 2K. T ef f( veI ) is not shown, 
as it is calculated from R vc \, which is normalised to -Rphot. 

Fig. 1181 shows the distribution of T e g amplitude for the 
three variations. The amplitude distributions of T e g^ pilot ) 
and T c ff( S p OC ) are both strongly peaked at lower values, but 
AT c ff( vcl ) shows a much wider distribution, and tends to- 
wards much larger amplitudes. The amplitudes of T e ff(phot) 
and T G ff( spGC ) mostly lie between 20 and 200K, but T e g( ve i) 
amplitudes tend to lie between 200 and 800K. The median 
amplitudes of the sample are 42.3, 72.56, and 582. 78K for 
and AT cff ( vcl ) respectively. Teff(vei), 
like T c ff( p hot), is a predicted temperature, whereas Teff(spec) 
is measured from our data. T e ff(vel) nas very large ampli- 
tudes in the context of radial pulsation. Adding this to the 
fact that T e ff(vei) does not agree with T e ff( spec ) suggests that 
radial pulsation can not explain the LSPs. 

Fig. [19] shows the relation of light amplitude to T e a 
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Figure 17. Distribution of median values of T ef j( phot ) (solid red 
line) and T e fj( spec ) (dashed blue line). 

amplitude for both T e ff( p hot) and T c ff( S p C c) ■ The graph shows 
that T c ft amplitude increases with increasing Mr amplitude, 
for both T eff{phot) and T eff(spec) . 



5 DISCUSSION AND CONCLUSIONS 

5.1 The Movement of the Visible Surface of a 
Sequence D Star during its LSP 

Possibly the most remarkable feature of the sequence D stars 
is the large change in the radial distance to the visible sur- 
face of a sequence D star as it passes through one LSP. By 
integrating the radial velocity with time we calculated in 
section 14.21 a radius change for radial pulsation of between 
3 and 6OR0, with the typical radius full amplitude being 
~41 R,0. The typical median radius of a sequence D star is 
around 135R.0. 

We can use this property to demonstrate just how un- 
likely radial pulsation is in these stars. An amplitude of 
41R.0 in a 135R© radially pulsating star corresponds to a 
fractional radius change of over 30 per cent from minimum 
to maximum radius. These radial amplitudes are very large 
in the context of radial pulsation: fundamental-mode Mira 
variables pulsate w ith radius changes of the order of 50Rq 
ijlreland et al.l [2004). but they have visible light amplitudes 
of ~6 magnitudes, while typical sequence D stars have light 
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Figure 18. Distribution of amplitude of effective temperature 
variation for T ef j( pllot ) (solid red line, forward shading), T ef j( spec ) 
(dashed blue line) and T e gy ve n (dot-dashed green line, back shad- 
ing). The higher end of the AT e fK ve n distribution has been cut 
off, in order to show greater detail at the peaks of AT^p^t) and 

ATeff(spcc) • 

amplitudes of AMr ^ 0.8 mag. It seems unlikely that such 
a modest light change could be associated with so large a 
radius change, suggesting that radial pulsation is not the 
cause of LSPs. 

This comparative analysis backs up a huge problem 
with the radial pulsation model which was raised in sec- 
tion [L2] the changes in radius would lead to changes in T e s 
that are vastly greater than the directly observed changes 
from spectra or photometric colour. 

In addition to the problems with radial pulsation raised 
by our own data, there are also other problems raised by 
IWood et all l|2004h : the length of the primary period does 
not vary during the LSP as expected; and to date there is 
no known mode of radial pulsation that has the required 
peri ods. 

IWood et al.l l|2004h proposed nonradial pulsation as a 
possible cause for the LSPs. They showed that the best 
match to the observed periods is given by the g mode with 
n — 2 and I = 2. In this case, the nett apparent radial veloc- 
ity arises from some parts of the visible surface approaching 
while other parts are receding. Thus at a given position on 
the stellar surface, the velocity amplitudes and overall radial 
motions will need to be even larger than in the radial pulsa- 
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Figure 19. Top Panel: AT^p^) plotted against AMr. Bottom 
Panel: AT e g( spec ) plotted against AMjj. 



tion case. Given the large amplitudes calculated in this work 
for radial pulsation, nonradial pulsation would lead to dis- 
tortions of the star that appear far too large to be credible, 
especially since g modes have signifi cant amplitudes only in 
radiative regions (|Wood et al.ll2004r ). 

It would therefore seem that the only way to generate 
the observed change in radial distance to the visible sur- 
face of a sequence D star is to have the star itself move, 
i.e. it must be in a binary system. Binary models, however, 
are inconsistent with various observational results described 
above. Most particularly, the non-uniform distribution of the 
angle of periastron and the low companion masses of these 
stars appear incompatible with a binary system model. The 
problem of the companion mass might be overcome if the 
companions were initially planets in eccentric orbits whose 
mass has been increased via accretion of material from the 
red giant. This would mean that although the companions 
currently have masses in the Brown Dwarf range, they were 
not born as Brown Dwarfs. To remove the problem of the 
non-uniform distribution of angle of periastron we would 
need to have an observational selection effect that rejects 
angles of periastron in the range 0°-180°, i.e. we would need 
an effect that masks the sequence D variation when perias- 
tron occurs with the smaller companion in front of the red 
giant. An effect which could produce this bias is difficult to 
imagine. 



As with past studies, we are unable to find definitive ev- 
idence for any model which explains the Long Secondary 
Periods in sequence D stars. To conclude, we therefore list 
all the currently known properties of LSPs. Any proposed 
model for LSPs should be able to explain all of these at- 
tributes. 

(i) Stars exhibiting LSPs occupy a clearly defined period- 
luminosity sequence. 

(ii) LSPs are of length -250-1400 days. 

(iii) LSP light variation is not regular and minima in par- 
ticular vary in depth from cycle to cycle (see Fig. 2 of Wood 
et al. 1999). 

(iv) The primary pulsation is visible in the light curve at 
all times throughout the LSP and the p rimary period doe s 
not significantly change with LSP phase l|Wood et al.ll2003 ) . 

(v) The ratio of LSP to pulsation period is ~8-10. The 
shorter period variation lies usually on sequence B and is 
thought to be the first or second overtone radial pu lsation 
(Wood et al. 1999; llta et ai]|2004 iFraser et al.ll2005l ). 

(vi) Lack of spectroscopic line broadening in observations 
of sequence D varia bles indicates that any rotational veloc- 
ities are «: 3km s" 1 l|Qlivier fc Woodll2003l ). 

(vii) The radial velocity curves show a characteristic 
shape: observed radial velocity increases quicker than it de- 
creases. 

(viii) The velocity amplitudes cluster tightly around 
3.5km s~ . 

(ix) The light-velocity phase shift for both the short- 
period and the LSP is ~0.25, with the light lagging the 
velocity. Minimum light is roughly aligned with mean rising 
velocity. 

(x) The ratio of colour to light variation of the L SP and 
prim ary oscillations are s imilar (IWood et alj 120041 ') . Simi- 
larly, iDerekas et al.l (|2006l ) showed that the ratio of blue to 
red amplitude of the LSP was similar to the ratio for stellar 
pulsation, and somewhat different from that due to ellip- 
soidal light variations. 

(xi) There is no correlation between LSP light and veloc- 
ity amplitude. 

(xii) Stellar radius tends to be between 100 and 2OOR0, 
with large radius variations of 3-6OR0 (if it is assumed that 
the LSP is caused by radial pulsation). 

(xiii) The equivalent width of the Ha absorption line 
varies with the LSP, indicating chromospheric activity 
i|Wood et al.ll2004| y 
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APPENDIX A: CALCULATING LUMINOSITY 
AND EFFECTIVE TEMPERATURE FROM 
PHOTOMETRY 

The luminosity was calculated from the MACHO Red pho- 
tometry. The MACHO Red and Blue magnitudes were first 
dered dened using a colour e xcess of E(B — V) = 0.08 for the 
LMC (|Keller fc Woodll2006l) and absorption A(\ ) calculated 
by the reddening law from ICardelli et al.l (|l989h . This gave 
A(M B ) = 0.2531 and A(M R ) = 0.1783. 

The bolometric correction to M Ro (where the zero sub- 
script denotes dereddened values) was derived through the 
following process. Note that we excluded the carbon stars in 
this procedure since the standard relations do not apply to 
them. First the bolometric correction to Vo wa s calculated, 
using a calibration derived from the data of iFluks et al.l 
l| 19941 ): 

BCv = 0.662 - 0.872(U - 7) - 0.128(U - 1)1, (Al) 

for 1 < (V — I)o < 4. The colour was converted from (Mb — 
Mr)q to (V — I)o using the equation 

nr n M Bo - M Ro - 0.120 
( ^- J) ° = 0625 

which was derived from MACHO magnitudes and V and 
I images of our field taken with the 40-inch telescope at 
Siding Spring Observatory, Australia. Substituting (V — I)o 
in equation [XT] gave BCv- 

The bolometric correction to Mb was calculated from 

BC Mb = BCv - {M B - V) 

where Vo was calculated using the formula 

Vo = M Bo + 0.042 - 0.155(M B - M R ) 

derived again from the 40-inch images. From BCm b it was 
a simple step to calculate BCm r using 

BCm r = BCm b — [Mr — Mb)o, 

the apparent bolometric magnitude was calculated from 
m bo i = Mr + BC Mr , 
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and the absolute bolometric magnitude was calculated from 
Mboi = m bo i - 18.54 

taking the distance modulus as 18.54 for the LMC (from 
iKeller fc Woodll2006h . 

Finally, the luminosity was calculated from 

L — Lq x lO '^ 4 - 76 -^"!) 

where 4.75 is the adopted absolute bolometric magnitude of 
the Sun. 

The effective temperature was calculated from the pho- 
tometry usuig_j 1 ijt_jri^eto_the stellar atmosphere derived 
data of iKucinskas et all (120051 ) for [Fe/H] = -0.5 (as ap- 
propriate for the LMC): 

3723 

T cff = for (V - I)o > 1.2. 

(( V-I)o- 0.8957)- 
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